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DIMENSION REDUCTION FOR DIPOLAR BOSE-EINSTEIN 
CONDENSATES IN THE STRONG INTERACTION REGIME 

WEIZHU BAO, LOIC LE TREUST, AND FLORIAN MEHATS 


Abstract. We study dimension reduction for the three-dimensional Gross- 
Pitaevskii equation with a long-range and anisotropic dipole-dipole interaction 
modeling dipolar Bose-Einstein condensation in a strong interaction regime. The 
cases of disk shaped condensates (confinement from dimension three to dimension 
two) and cigar shaped condensates (confinement to dimension one) are analyzed. 
In both cases, the analysis combines averaging tools and semiclassical techniques. 
Asymptotic models are derived, with rates of convergence in terms of two small 
dimensionless parameters characterizing the strength of the confinement and the 
strength of the interaction between atoms. 


1. Introduction and main results 


In this paper, we study dimension reduction for the three-dimensional Gross- 
Pitaevskii equation (GPE) with a long-range and anisotropic dipole-dipole interac¬ 
tion (DDI) modeling dipolar Bose-Einstein condensation [11, 14]. In contrast with 
the existing literature on this topic [ 1 ], we will not assume that the degenerate 
dipolar quantum gas is in a weak interaction regime. 

Based on the mean field approximation [3, 9, 13, 18, 19, 20], the dipolar Bose- 
Einstein condensate is modeled by its wavefunction T := ^(f, x) satisfying the GPE 
with a DDI written in physical variables as 


ihdt* = + U(x)tf + Ng\*\ 2 * + AC dip ([/dip * I'Ll 2 ) , (1.1) 


where A is the Laplace operator, V(x) denotes the trapping harmonic potential, 
m > 0 is the mass, h is the Planck constant, g = i7Th rn as describes the contact (local) 
interaction between atoms in the condensate with the s-wave scattering length a s , 
N denotes the number of atoms in the condensate, and the dipole-dipole interaction 
kernel [/d; p (x) is given as 


1 l-3(x-n) 2 /|x| 2 

u *> (x) - a -b#- 


1 1 — 3cos 2 ($) 
4-7T Ixl 3 


x e 


( 1 . 2 ) 


with the dipolar axis n = (ni,ri 2 ,n 3 ) € M 3 satisfying |n| = \/n\ + n\ + n 3 = 1 . 
Here 6 is the angle between the polarization axis n and relative position of two 
atoms (that is, cos 9 = n • x/|x|). For magnetic dipoles we have Cdi p = /D/^dip’ 
where go is the magnetic vacuum permeability and //dip the dipole moment, and for 
electric dipoles we have Cdi p = Pdi p / e 0i where eo is vacuum permittivity and pdi p 
the electric dipole moment. The wave function is normalized according to 


| (t, x)| 2 dx = 1. 


1 






2 


W. BAO, L. LE TREUST, AND F. MEHATS 


1.1. Nondimensionalization. We assume that the harmonic potential is strongly 
anisotropic and confines particles from dimension 3 to dimension 3 — d. We shall 
denote x = (x,z), where x G R 3 ~ d denotes the variable in the confined direction(s) 
and z G R^ denotes the variable in the transversal direction(s). In applications, 
we will have either d = 1 for disk-shaped condensates, or d = 2 for cigar-shaped 
condensates. Similarly, we denote n = (n x ,n z ) with n x G R 3 ~ d and n z G R d . The 
harmonic potential reads [2, 15, 16] 

n«o = y(^M 2 +^m 2 ) 

where oj z S> uj x . We introduce three dimensionless parameters 


— V Wi/ ^zi P — 


4irN\a s 


Aq — 




dip 


ao 3|g| ’ 

where the harmonic oscillator length is defined by [2, 15, 16] 

°° \mw x ) 

The dimensionless parameter Ao measures the relative strength of dipolar and s- 
wave interactions. Let us rewrite the GPE (1.1) in dimensionless form. For that, 
we introduce the new variables t, x, z and the associated unknown T defined by 

t = u x t, x= — , 5 = —, 'ffit, x, z) = aJ f (f,x, z). (1.3) 

a 0 a 0 

The dimensionless GPE equation reads [2, 15, 16] 


idfi L = —-A'k + - 
* 2 2 


I~|2 i 1 i ~i2 

N + ? M 


T + /3ct|T| 2 T + 3A 0 /3 U d ip * |4| 2 *, (1.4) 


where a = sign a s G {—1,1}. Define the differential operators d n = n • V and 
<9 nn = d n d n . Mathematically speaking, the convolution with U d j p in equation (1.1) 
has to be considered in the distributional sense and we have the following identity 
(see [3]) 


Ldip(x) = p.v. 


1 


47T | X |' 


1 - 


3(x • n) 


x 


= -^(x) -d n 


1 


47r|x| 


x G 


(1.5) 


with S being the Dirac distribution. 

Remark 1.1. Let us define the Fourier transform of a function u € L 1 (M 3 ) by 
u(k) = f u(x)e _ * fc ' x dx, x € R 3 . 

Jr a 


From identity (1.5), we get 


U d ip{k) = ~ ^ , for all k € R 3 . 


( 1 . 6 ) 


We can re-formulate the GPE (1.4) as the following Gross-Pitaevskii-Poisson 
system (GPPS) [3, 7] 

id$ = —AT + \ (jz| 2 + ^rl 5 ! 2 ) * + - Ao)|$| 2 $ - 3A 0 /?(d nn ^) 4, 

Aip = — |T| 2 , lim </?(f,x) = 0. 

|x|—>-oo 


(1.7) 
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Under scaling (1.3), dimension reduction of the above GPPS (1.4) was formally 
derived from 3D to 2D and ID in [1, 7, 17] for any fixed f3, Ao and n when e —>• 0 + . 
Rigorous mathematical justification was only given in the weak interaction regime, 
i.e. when (3 = 0(e) from 3D to 2D and when (3 = 0(e 2 ) from 3D to ID [1], It is an 
open problem for the case where (3 is fixed when e —>• 0 + . 


1.2. New scaling. In order to observe the condensate at the correct space scales, 
we will now proceed to a rescaling in x and Let us denote 

a = £ 2d / n f3- 2 In. (1.8) 

The scaling assumptions are 

a<l and e<Cl, 

We define the new variables 


t = t, 


z' = 


x' = a 1 // 2 x, 


which means that the typical length scales of the dimensionless variables are e in the 
z-direction and a -1 / 2 in the ^-direction. The wavefunction is rescaled as follows: 

y £ ’ a (t',x',z') := £ d/2 a- n/A ^(t,x,z)e iid/2£2 . 

Notice that the L 2 norm of T £,a is left invariant by this rescaling, so we still have 


/ 


\^ £ ' a {t,x,z)\ 2 dxdz = 1 . 


We end up with the following rescaled GPE (for simplicity we omit the primes 
on the variables): 


qj nr' \x\ 2 

iad t ^ £ ' a =—n z ^ £ ’ a - — A X T £ ’“ + 




+ a[a\V £ ’ a \ 2 + 3\ 0 pU^*\V £ 

where the transversal Hamiltonian is 


e,«l 2 j xJ/ £ > Q 


(1.9) 


1 . |z | 2 d 

U z := —A, + —- 

2 2 2 


and Ufa “ is defined by 


U e d g(x,z) = U& 


ip 


a 


,£z) , (x,z)G 


Let us remark that 


uS(k x , k z ) = £~ d a n / 2 Ud\p ( — ) , for all (k x , k z ) € M 3 


( 1 . 10 ) 


Thanks to identity (1.6), we can remark that U &; p is a bounded function of R ' 3 into 
[—|,|], For 7 > 0, we denote by V2 p the tempered distribution whose Fourier 
transform is 

1 . ( 7 k x ■n x + k z ■ n z f 


Vg p (k x ,k z )=(-- + 


\^k x \ 2 + |fc 2 | 2 


so that V dip (k x ,k z ) € [—1/3, 2/3] for all (k x ,k z ) € R 3 and 

&(k x ,k z ) = £~ d a n / 2 v£ e (k x ,k z ), for all (k x ,k z ) € R 3 . 


( 1 . 11 ) 
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Let us note that (1.8) is equivalent to 

P£~ d a n / 2 = 1 

so that equation (1.9) becomes 


2 | |2 

iad t ^ £ ' a =%U z ^ e ' a - — A a .T £ ’ a + L^T £ ’ Q 


+ a ( u |^ £ ’“| 2 + 3 \ 0 V ^ £ * | T £ ’ a | 2 ) ^ £ ’ a . ( 1 . 12 ) 

Remark 1.2. The spectrum ofTL z is the set of integers N. We define (oJk)k£N an 
orthonormal basis of L 2 (M 3 ) made of eigenvectors of 7~L Z where ujq is the ground 
state (associated to the eigenvalue 0) 

u 0 (z) = 7 ,r- rf / 4 e -l*l 2 / 2 . 

Remark 1.3. Since (L^p^x) is uniformly bounded in L°° and 

Lip -+ iS P «.e. 

as 7 —>• 0, Lebesgue’s dominated convergence Theorem ensures that 

V^* U ^V° P *U in L 2 (M 3 ) 
for all U € L 2 (M 3 ). Moreover, let us remark that 


3 d 


V^*U(x,z) = 


U(x, z ), (x, z ) € 


for all U such that U(x,z ) = V(x, \z\) for all (x,z) € 


In this paper, we study the behavior of the solution of equation (1.12) as e -A 0 
and a -A 0 independently so that /? may be bounded but can also tends to +oo. 

Our key mathematical assumption will be that the wavefunction 'k e,a at time 
t = 0 is under the WKB form: 

^ £ ’ a (0, x, z) = ^f nit (x, z) := Aq(x, z)e iS °W a , V(®, z) € M 3 . (1.13) 

Here Aq is a complex-valued function and Sq is real-valued. 

Let us introduce another parameter 7 > 0 to get a better understanding of the 
different phenomena involved during the limiting procedures. In this paper, we will 
study instead of equation ( 1 . 12 ) the following one : 

iadtfi = fiTL z i( - ^-A x if + kjrfi + a (a\fi\ 2 + 3A 0 V2 p * \fi\ 2 ^j if, (1.14) 

if(0,x,z) = Aq(x, z)e lS °^ a , \/(x,z) € M 3 . 

From now on, we denote by \E ' £,a ’ 7 the solution if of equation (1.14). Let us insist 
on the fact that \H £,a ’ 7 equal to the solution 'F £, “ of equation ( 1 . 12 ) if we assume 
that 7 = £y/a. 

1.3. Heuristics. In this section, we derive formally the limiting behavior of the 
solution of (1.14) as e (strong confinement limit), a (semiclassical limit) and 7 
(limit of the dipole-dipole interaction term) go to 0. Our main result, stated in the 
next section, will be that in fact these limits commute together: the limit is valid 
as s, a and 7 converge independently to zero. Thus, this gives us as a by-product 
the behavior of the solution of equation ( 1 . 12 ) as e and a converge independently 
to zero. 
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a) Strong confinement limit : e —>• 0. Let us fix a € (0,1] and 7 € [0,1]. Following 
[6], in order to analyze the strong partial confinement limit, it is convenient to 
begin by filtering out the fast oscillations at scale e 2 induced by the transversal 
Hamiltonian. To this aim, we introduce the new unknown 

$ £ ’ Q ’ 7 (t, •) = e itn z /e 2 ye,a,-y^ 

It satisfies the equation 

nfi \ x \ 2 / f 

iadt $ £ ’“’ 7 = A . c $ £ ’“’ 7 + ^$ £ ’°’ 7 + aF 'y —, $ £ ’ Q ’ 7 

2 2 V^ 2 

where the nonlinear function is defined by 

F 7 (0,$) = e l()H > (v\e~ im ^\ 2 + 3\ 0 V2 ip * \e~ im ^\ 2 ^ (1.15) 

A fundamental remark is that for all fixed <f>, the function 6 1 —)■ i ?7 (0,< h) is 2n- 
periodic, since the spectrum of 7i z only contains integers. For any fixed a > 0 and 
Ao = 0, Ben Abdallah et al. [6, 5] proved by an averaging argument that we have 
<!>£:«,7 = <|> 0 ,a ,7 _|_ where 4> 0,a ” 7 solves the averaged equation 

iad t <S>°’ a ^ = ~ A X $°’ Q ’ 7 + l|1$°-“- 7 + a F2 v ($°’ Q ’ 7 ), 4>°’"’ 7 (f = 0) = T 


a 

init 


(1.16) 

where F2 V is the averaged vector field 

F 2 v m = ^J^ F1 {.e^)de. (i.i7) 

In our study, we consider the case Ao € M and a similar averaging argument should 
give us the same result 

4 > £ ’ Q ’ 7 = $ 0 ,a ,7 + 0(s 2 ). 

b) Semi-classical limit : a —>• 0. Let us remark that equation (1.14) is written in the 
semi-classical regime of "weakly nonlinear geometric optics", which can be studied 
by a WKB analysis. Here we are only interested in the limiting model, so in the first 
stage of the WKB expansion. Let us introduce the solution S(t,x ) of the eikonal 
equation 


5*5 + 


|V a S|“ \x 


+ — = 0, 5(0, x) = S 0 (x) 


and filter out the oscillatory phase of the wavefunction by setting 

^E,a,7 _ g -iS(t,x)/a ^£,£*,7^ 

so that 


1 


d t n £ ’ a ^ + v x s ■ v x n £ ’ a ’^ + ^n £ ’ a ^A x s = A x n £ ^ - i^-n £ ’ a ^ 

-i (u|H £ ’ Q ’ 7 | 2 + 3A 0 F d ]p * |H £ ’ Q ’ 7 | 2 ) H £,a ’ 7 , 

where 

O e,a, 7 ( 0 , x, z) = Ao(x, z), for all (x, z) € M 3 . 

For all fixed e > 0, we can expect that 

H £,a ’ 7 = H £ ’ 0,7 + O(a), 


,n z 


(1.18) 

(1.19) 

( 1 . 20 ) 
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as a —>• 0 where fl e ’ 0,7 solves the equation 

d t n £ ’°’ 7 + X7 x S • V x ft £ ’ 0 ’ 7 + yfl £ ’°’ 7 A X S (1.21) 

= -i^O e ’ 0 ' 7 - i (a|ft £ ’°’ 7 | 2 + 3A 0 F d ] p * |fl £ ’°’ 7 | 2 ) fl £ ’ 0 ’ 7 , 

n £ ’°’ 7 (0, x, z) = Ao(x,z), for all (x,z) G M 3 . 

Remark 1.4. A key point here in this analysis is that the nonlinearities l 77 and 
Fav are gauge invariant i.e. for all U G F 2 (R 3 ), all 7 G [0,1] and for all t, we have 

F 7 (f, Ue iS ' a ) = F 7 (t, U)e iS / a , F 7 v (Ue iS / a ) = F 7 v (U)e iS / a . 

c) Dipole-dipole interaction limit 7 —>• 0. We expect that for any (e,a) G (0, l ] 2 

_ vj/e,a,0 + Ofa q ) 

where q > 0 and 


2 I 12 

iad t y £ ’ a ’° = ^n z ^ £ ’ a ’° - yA ^^’ 0 + 

+a (<r|tf e -“' 0 | 2 + 3AoF d ° ip * |^ £ ’"’°| 2 ) ^ £ ’ Q ’°, 
* £ ’ a ’°(t = 0 ) = ^ nit . 


( 1 . 22 ) 


In this paper, the main difficulty we have to tackle and also the main difference 
with respect to the previous work of the authors [4] in the case Ao = 0, is the study 
of this limit 7 —>• 0 . 

d) The simultaneous study of the three limits. We introduce for any (e, a, 7 ) G (0, l ] 3 
A £ ’“’ 7 (f, X , z) = F^/e\-iS{t,x)/a ^ ^ for ( x> z ) € M 3 , 

which is the solution of the equation 

<M £ ’ a ’7 + V x 5 • V X A £ ’ Q ’ 7 + ^A £ ’ a ’ 7 A x S = ^^A £ ’ a ’ 7 - iF 7 ^-L,A £ ’ a ’^ ,(1.23) 

A £,a ’ 7 (0, x, z) = Ao(x, z). 

We will also consider the solution A £ ’ 0,7 of (1.23) with a = 0, the solution A £,a, ° of 
(1.23) with 7 = 0 and the solution A 0, “’ 7 of 

8 t A °’ a ’ 7 + V X S ■ V x A°’ a ’ 7 + yAL°’“’ 7 A ;c = ^^A °’“’ 7 - iF 7 v (A°’ Q ’ 7 ) , (1.24) 

A°’ Q,7 ( 0,x,z) = Aq(x,z), 

for all (x,z) G R 3 . As long as the phase S(t , •) remains smooth, i.e. before the 
formation of caustics in the eikonal equation (1.18), we expect to have 

A £ ’ a ’7 = A 0 ’ 0 ’ 0 + 0(e 2 + a + 7 9 ), 

and the solution d /£,a ’ 7 of equation (1.14) is expected to behave as 

4r e ’ a ’ 7 (t, x, z) = e -m z /e 2 e iS(t,x)/ a A 0,0,0fa Xj z ) + 0(e 2 + a + 7 9 ) (1.25) 

for some q > 0 . 
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1.4. Main results. In this paper, our main contribution is the rigorous study of the 
dipole-dipole interaction limits 7 —>• 0 as well as the study of the three simultaneous 
limits e -A 0, Oi —^ 0 and 7 -A 0 involved in the problem. The techniques used for 
the study of the limits e —>• 0 and a —> 0 were developed by the authors in [4]. We 
will recall and use some of the results proved in this first paper. 

1.4.1. Existence, uniqueness and uniform boundedness results. Let us make precise 
our functional framework. For wavefunctions, we will use the scale of Sobolev spaces 
adapted to quantum harmonic oscillators: 

B m (R 3 ) := {u € H m (R 3 ) such that (| x\ m + \z\ m ) u £ L 2 (R 3 )} 

for m £ N. 

Remark 1.5. Assuming that m > 2, we get that 

B m (R 3 ) -A H m (R 3 ) -A L°°(lR 3 ). 
are two algebras. In this paper, we will also make frequent 


(1.26) 


H m (R 3 ) and B 
use of the estimate 

|||x| fc < 9>|| L 2 < C'HuHsm, for all u £ B m (R 3 ) and k + \n\ < m 
(see [12] and [6] for a more general class of confining potential). 

For the phase S, we will use the space of subquadratic functions, defined by 
SQ k (R 3 ~ d ) = {/ £ C fc (R 3 " d ;R)) such that dff £ L°°(R 3 " d ), for all 2 < |k| < k}. 

where k £ N, k > 2. In the following theorem, we give existence and uniqueness 
results for equations (1.18), (1.23) and (1.24), as well as uniform bounds on the 
solutions. 

Theorem 1.6. Let (e,a, 7) € [0,1] 3 , Ao £ B m (R 3 ) and Sq € SQ s+ i(R 3-rf ), where 
m > 5 and s > m + 2. Then the following holds: 

(i) There exists T > 0 such that the eikonal equation (1.18) admits a unique 
solution S € C([0,T];SQ s (R 3 - d )) nC s ([0,T] x R 3 ~ d ). 

(ii) There exists T £ (0, T] independent of e, a and'y such that the solutions A e,a,7 
and A 0 ' 11 ' 1 of, respectively, (1.23) and (1.24), are uniquely defined in the space 

C([0,T] ; r(R 3 ))nC 1 ([0,T];fi m - 2 (R 3 )). 

(in) The functions (A £ ' a ' 1 ) e , ct ^ are bounded in 


3 )) 


C([0, T]; B m {R 3 )) n C 3 ([0, T]; B m ~ 2 ( 

uniformly with respect to (s,a, 7 ) £ [ 0 , l] 3 . 

1.4.2. Study of the limits a —» 0, e —> 0 and 7 —>• 0. We are now able to study the 
behavior of A e,a,7 as a —>■ 0, e —>• 0 and 7 ^ 0 . 

Theorem 1.7. Assume the hypothesis of Theorem 1.6 true. Then, for all (e, a, 7 ) £ 
[ 0 , l] 3 , for all q £ ( 0 , 1 ), we have the following bounds: 

(i) Averaging result: 


|Me,a ,7 _ ^40,0,71 

Semi-classical result: 

WA 6 ’ 01 ' 1 — y 4 £ ’°’'i| 


L°°([ 0 ,T];S m - 2 (R 3 )) — 


< Ce 2 


L°°([ 0 ,T];S m - 2 (R 3 )) — 


< Ca 


(1.27) 

(1.28) 
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(Hi) Dipole-dipole interaction limit result: 

ll^- £ ’ a ’ 7 — ^ £ ’ a ’°llL°°([0,T];_B m - 5 (K 3 )) — (1-29) 

(iv) Global result: 


I AS,a,y _ 40,0,0m _ c o < /7 ("=-2 , , 

I 71 71 HL°°([0,T];B m - 5 (R 3 )) — + a 7 ). 


(1.30) 


The constants C and C q do not depend on a, £ and 7 but C q does depend on q. The 
estimates related to the original equation ( 1 . 12 ) can be summarized in the following 
diagram: 



Remark 1.8. The case Ao = 0 has already been studied by the authors in [4] where 
we got estimates that are similar to (1.27) and (1.28). 

Remark 1.9. Assume that either n x = 0 or n z = 0. Then, for all (s,a) € (0, l] 2 , 
for all q such that 

{ 9 = 1 ifd = 1 , 

1 <? e [1,2) i/d = 2, 

we get the same conclusion as in Theorem 1.7. 

The following immediate corollary gives a more accurate approximation of 
j^£,a,£y/a than A 0,0,0 . This result can be useful for numerical simulations and has to 
be related to the ones of Ben Abdallah et al. [5]. 

Corollary 1.10. Assume the hypothesis of Theorem 1.6 true. Then, for all (e, a) € 
[ 0 , 1 ] 2 ; we have the following bound: 

< C(e 2 +a) 

where C > 0 does not depend on e or a. 

The following proposition concerns the special case of an initial data polarized 
on one mode of Ti z . It generalizes the case studied by Bao, Ben Abdallah and Cai 
[1, Theorems 5.1 and 5.5] where the initial data was taken on the ground state of 
U z . 

Proposition 1.11. Let k € N. Assume the hypothesis of Theorem 1.6 true. Assume 
also that 

A 0 (x, z ) = a 0 (x)uk(z), (x, z) € M 3 

where is defined in Remark 1.2. Then, the function A 0,a ’ 7 stays polarized on the 
mode uih i.e. 

A°’“ , 7 (t, x, z) = B a ^(t,x)u k (z) for all z € R d . 
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Here, F Q ’ 7 is the solution of 

d t B a ’ 7 + V X S ■ V X F Q ’ 7 + i B a ’" / A x S = ^^F "’ 7 - fG 7 ’ fe (F Q ’ 7 ), (1.31) 


where 


Let 


F“’ 7 (0, x) = ao(x), x 6 M 3 

G 7 ,fc (u)(x) = u(x) J \u k (z)\ k (a + 3A 0 F d 7 p * \u) k u\ 2 (x, z)^j dz. 


j q = 1 if d = 1 

\qe[l,2) ifd = 2, 

we have moreover the following bound for all a £ [0,1] : 

P 0,a,7 40,a,0ii _ ^ n 

A ll/ oo ([0,T];B m — 5 (R 3 )) — 

where C q does not depend on a but depends on q. Hence, we obtain that 
|| jse,a,ey/a _ J 4 0 > 0 >°|| ioo ^ () _ . Bm -s^3^ < C(s 2 + a + (£y/a) Q ) 
and for a £ (0,1] fixed 


\q,£,a,£y/^ _ ^ 


O.o.O | 


< Ge«. 


(1.32) 


lL°°([0,T];B m - 5 ( 

Remark 1.12. Let us notice that by Remark 1.3, the nonlinearity G 7 ’^ of equation 
(1.31) becomes a local cubic nonlinearity when 7 = 0 

" n 2 — d, 


G°’ k (u)(x) = ( Uz 3d d \u(x)\ 2 u(x), for all x £ 


p3 —d 


The paper is organized as follows. In Section 2, we study some properties of the 
dipolar term that are needed in the proofs of Theorems 1.6, 1.7 and 1.11 given in 
Section 3. 


2. Study of the dipolar term 

Let us define for 6 £ M, 7 € [0,1] and <h £ L 2 (R 3 ) 

Fi(9, <L) = e im * (a\e- im ^\ 2 e- im ^) , 

F ltao {*) = ±fi* Fl (6,Q)M, 

F 7 ( 0 ,$) = 3\ 0 e ieH ° (V 7 p * |e-^$| 2 ) e - i6H *<&, 

& J? 

so that 

F 1 = F\ + F 7 and F 7 = F hav + F 7 ^. 

I 11 order to prove the uniform well-posedness of the nonlinear equations (1.23) 
and (1.24), we will need Lipschitz estimates for F 7 (0, •) defined by (1.15) and F,7„(-) 
defined by (1.17). We only study here the dipolar terms F^ and Fff av since the 
cubic ones Fi (#,<!>) and Fi tav ($>) have already been studied in [4, Lemma 2.7.] (see 
also [ 6 , Proposition 2.5], [10, Lemma 4.10.2] or [ 8 , Lemma 1.24]). 
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2.1. Some properties of K/(#, •) and Using the fact that takes its 

values in [—1/3, 2/3] (Remark 1.1), we get the following lemma. 

Lemma 2.1. Introduce the convolution operator 

K 1 : u € H m ( M 3 ) i—► R d ] p *u£ H m (R 3 ) 

for 7 € [0,1] and m G N where V[/ p is defined by (1.11). Then, we get for all 
u € LT m (M 3 ) 

The following lemma gives Lipschitz estimates for the dipolar terms. 

Lemma 2.2. For all m > 2 and M > 0, there exists C > 0 such that 
\Wlaviu) - Fl av {v )|| Bm < UM 2 ^ - 
||^ 2 7 (<?,«) - ^ (MII B m < UM 2 ||r - W || Bm , 
for all u, v € B m (R 3 ) satisfying ||r||_b™, < M, ||u|| Bm < M, for all 8 G R and for all 

7 € [ 0 , lj. 

Proof. Let us fix 7 € [0,1], u, v € B m (R 3 ) satisfying ||u||_Bm < M, ||u||sm < M. To 
begin, assume that 6 = 0, then we get that 

117 (0,u) - F 7 (0,b)|| b „ = 3A„ \\K\\u\ 2 )u - K 7 (|t>|>|l B .» 

< 3Ao ||A 7 (M 2 )(» - „)|| Bm + 3A„ ||(A 1 '(|ti| 2 ) - 1C{ M 2 )H| b „ . 
Lemma 2.1 and Remark 1.5 ensure that 


\\K n (\u\ 2 ){u-v)\\ Hm < C\\K^(\u\ 2 )\\ Hm \\u-v\\ H m <C\\[u\ 2 \ 

< C \\u\\ 2 H m \\u — v\\ H m < CM 2 \\u — v\\ B r, 


u — v\ 


We also have 


x 


‘ Jr 7 (l“l 2 )(«-®)|| i i < C||K 7 (M 2 )|| b „ |||x|"*(u-»)|| i 2 


< C ||A' 7 (|ll| 2 )|| Bm ||u - u|| B m < CM 2 ||u - ti|| B „ 
For the second term, we get 


||(^(|u| 2 ) - K^\v\ 2 ))u\\ Hm < C\\K^\u\ 2 - \v 


H m II^IIh” 


<c||h 2 - 


v\ \\ H m IMItfm < CM 2 ||U — v\\ B n 


and 


x 


m (KU\u\ 2 ) - KU\v\ 2 ))u\\ l , < C |j A 7 (|u | 2 - M 2 )|| b „ |||xr«|| B , 
<C||H 2 - 


U | || fjm 11 ^ 11 B rn — 11 U ^ 11 


This gives us 


\\Fl (0,«) - Fl (0, u )|| Brn < CM 2 ||u - u| 


B r, 


where C depends on m but is independent of 7 , u and v. Since are isometries 

of B m , we get for 6 £ M 


\\F2(e,u)-F2(9,v)\\ Bm = 


iOHz ( zp'y ( n o—iSH 


7 I n o—iSH 


< 


FI (0,e- iOH 'u) -Fno,e 


7 in c—iOtiz 


Ff [0,e- wn ~~u) - Ff ( 0, e 


B" 


< CM 2 


2 


-iOHz 


V 

(u — v) 


B n 


B” 


< CM 2 llu-d 


B" 
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and 


= ^J a (*?(»,<*)-*?(«,- v))d« 
< ^ J* \\I? W,u) - F? ($, 


Bm 


v)\\ B m dd 


< CM 2 Ilit — v\ 


B m 


□ 


2.2. The limit 7 —>• 0. Let us study now the behavior of and av as 7 —>• 0. 
2.2.1. General case. 

Lemma 2.3. For all m > 2, q € (0,1), there is a constant C m _ q > 0 independent 
of 7 such that 

||F 7 (9,u) - F° (6,u) \\ Bm < 7 9 C'm,g||w|||m||u|| B m+5, 

II F av ( u ) - F av (“)|| B m ^ l q Cm,q\\u\\ 2 B m\\u\\ B m+ 5 , 

for all u € I? m+ 5 (M 3 ), for all 7 € (0,1] and for all 9 € R. 

Proof. Let u € R m+ 5 (R 3 ) and 7 € (0,1]. As in the proof of Lemma 2.2, we can 
assume that 9 = 0. Thanks to Remark 1.5, we get 

||F 7 (0, u) - F° (0, u)|| Bra = ||in(M> - K°(\u\ 2 )u\\ Brn 

<c Y \\d K (K\\u\ 2 )u- K\\u\ 2 )u)W l2 

\n\<m 

+C|||xr (IO(\u\ 2 )u- K°(\u\ 2 )u)W l , 

<c Yl ll- ff7 ( 8 '"M 2 ) - ^ 0 (8 m I“I 2 )|| l » lld'°“lli’ 

|ki|+|«s| <m 

+C\\K\\u\ 2 )-K\\u\ 2 )W 1 

<C\\u\\ B m Y I (^dip - ^dip) *(d K \ u \ 2 ) 

|«| <m 

Let us denote v := d K \u\ 2 for some |k| < m. We have that v € R 5 (R 3 ) and, since 

B m ^ L °°, 

||u|| B 5 F 11 U U^Jm+5 F G 11 U 11 11 U 11 B m+5 . 

Moreover, 


\-x\ m u\\ L 2 




< C 


L°° 


Ulp - LSp 1 v 


Li 


Since for all (k x , k z ) € 


^dip(kx> k z ) — i o + 


1 . ( 7 k x ■n x + k z ■ n z f 


\"yk x \ 2 + \k z \ 2 


Vg p - ^di P = W?(k x , k z ) + W2(k Xl k z ) 


we obtain 
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where 


W?(k x ,k z ) = - 
W 2 {k x ,k z ) = 


(k z ■ n z ) 2 \'yk x \ 2 7 2 {k x ■ n x ) 2 


i\lk x \ 2 + \k z \ 2 ) \k z \ 2 It^xI 2 + \k z \ 2 
27 (k x ■ n x )(k z ■ n z ) 


Step 1 : Study of 


W?v 


L 1 


hk x \ 2 + \k z \ 2 
. Since we have 


\W?(k x ,k z )\< l7fca;|2 


(hk x \ 2 + \k z \ 2 ) 1 + J^ 4 ’ 


we get for q\ > 0 and p\-/p\ £ [ 1 , +00] such that = 1 , by Holder, 


\W^ 


L 1 


< 


f 

J R E 


1 + 


I k z 


-dk x dk z < ||/i|| L p' IIsiIIlpi 


\lk x 


where 


fi = h\k x \) d/pi (l + \k x \Y \v 

1 

9 i = - 


{ , y\k x \) d / pi (1 + \k x \ 2 ) qi ( 1 + 


I k z 


|7 k x 

Thanks to the change of variable k = k z /\'yk x \, we obtain that ||<?i||lpi does not 
depend on 7 : 

H npi _ f dk x f dk 

ll9lllLP1 L-«(l + \k x \ 2 ) piqi L (l + |fc|T 

so that \\g\ \\ r ^ Pl < +00 if and only if 

PiQi > (3 — d)/2 and pi > d/2. 

On the other hand, we have 

ll/i|l ir ; = \\\Kf r ' (1 + IfelT ISIII^;. 

Assume that pi € [1,2) n (d/ 2,2) is fixed and define m, — 3 ( y j. Then, we get 
thanks to Sobolev inequalities and inequality (1.26) of Remark 1.5 that 

ll/illj/, <cv </p ‘llfel' f/ '’ 1 0 + IMT <c 7 J /'"||f|l B „ 1+ aL«„. 

In the case d = 1, we choose p = pi = 1 and q\ £ (1, 5/4) so that 

d 

mi H-b 2gi < 5. 

Pi 

In the case d = 2, we fix p £ [1,2). Then, we choose pi = 2/p £ (1,2] and 
(/i £ (l/2pi, 3/4) so that mi £ [0, 3/2) and mi + /y + 2</i < 5. We proved that 


for 


W/V < 7 P Cp||u||Bm ||u|| B m+5 


[ p = 1 if d = 1 

[ p £ [1,2) if d = 2 . 
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Step 2 : Study of 


Wjv 


L 1 


We have 


\Wj(k x ,k z )\ < 


, IM 

J l\k x \ 


1 + 


%\ 2 

l 2 \k x \ 2 


so that, we get for 02 > 0 and P 2 ■ Vo € fl, + 00 1 such that-b - 7 - = 1, 

P2 P 2 


\w2(kx,k z )\\v\ 


L 1 


< 


i 

J RS 


2\v\ 


\kz\ 

l\k x \ 


1 + 


1^1" 


■dk x dk z < 2\\f 2 \\ LP ' 2 \\g2\\LP2 


7 2 I*x| 2 


where 


f2 = {l\k x \) d/P *(l + \k x \ 2 ) q2 \V\ 

1*2 I 


92 = 


l\kx\ 


{n\k x \) d/P2 (1 + |^| 2 ) 92 (1 + 


l*2| 2 

|7**I : 

Thanks to the change of variable k = p^r-p we get that | g 2 \ \ rj ‘2 does not depend on 

T 

11 ups = f dk x f \k\ P2 dk 

imiLP2 J R3 - d (i + \k x \r 292 L {i + \k\T 

and \\g2\\ P LV2 < +00 if and only if 

P 2 Q 2 > (3 - d)/2 and p 2 > d. 

Let us fix from now on q € (0,1) and P 2 = d/q € (d, + 00 ). For p% > 0 and <73 > 0, 
we get thanks to Remark 1.5 that 


Wf2\\ L P 2 =i q 

< 'yl 

< C 7 9 


feb /w (i + fei 2 ) te+,B) (i + iLiTi«i 


L p 2 


(i + \k x \ 2 ) q3 (i + \k z ^r 


(i+ifei j )®(i+ifc.i 2 r 

1 


L p 2 


P 3 


and 


(1 + |**| 2 )® , (1 + |*2| 2 ) M 

1 


, IMI_B3/2+d/p2+2(<J2+93+P3) 


L p 2 


L p 2 


(i + |**| 2 ) 9S (i + |*2| 2 ) P3 

if and only if 

293 P 2 > 3 — d and 2psp' 2 > d. 

Let us choose 

/ (3 - d) (3-d) 1 \ 


< 00 


(3 ~ d)(p2 — 1) (3 — d){p2 — 1) 1 
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and 


so that 


P3 G 


( d(p2 - 1 ) d(p 2 - 1 ) J_ 
V 2p 2 ’ 2p 2 + 12 


3/2 + d/p 2 +2(q 2 + q 3 +p 3 ) < - + - + 
and 


3 1 d + (3 - d) + (3 — d)(p 2 - 1) + d(p 2 — 1) 


P 2 


< 5 , 




L 1 


< l q C q \\u\\ gm || u\\ £?m+5 . 


Hence, 

||F 7 (0,7/) - F° (0,7i)|| sm < 7 9 C' mi9 ||7/||| m ||7/|| B m+5. 

Finally, we get that 

\\ F 2 v («) - («)|| Bm < h lo ll^ 7 (M) - F° (M)|U de < 7 , C , m) ,||«||| tn ||u|| fl m+6. 

□ 

2 . 2 . 2 . Case of a function which is polarized on one mode ofTi z . 

Lemma 2.4. Let m > 2, M > 0 and 

q = 1 if d = 1 , 

q e [1,2) 7 /d = 2. 

Then 

\\ F av («) - ^ («)|| B m < 7' ? C' mi g||7/||| m ||7/|| S m+5, 

/or all 7 € ( 0 , 1 ] and all u € H m+ 5 (M 3 ) under the form 

u(x,z ) = a(x)ujk{z), for all (x,z) € M 3 

where k € N and uj k is defined in Remark 1.2. The constant C m , q depends neither 
on u nor on 7 . 

Proof. Let 7 € [0,1] and u €E H m+ 5 (R 3 ) such that u(x, z) = a{x)ujk{z). We get that 

F 2v( U )= U} k / UJ k (z)F 1 (0,u){-,z)dz +U k / (jJk{z){F 2 l + F 1 2 ) ' 2 )(0, t/)(-, z)dz 
J K d .7R d 

where 


F 7 ,(M) = 3A 0 e ieWz f W 7 * \e~ im *u\ 2 ) e~ tm *u 


-i 8 H z „.\ 2 \ —idP. 


F 2 7 2 (d, u ) = 3A 0 e*^ ( W 7 * |e 




u ~ e 


( 2 . 2 ) 

(2.3) 


for all u € L 2 (R 3 ), 7 € [0,1] and # € R. We also have 0 Jk(z) = ±Wfc(—z) for all 
z € R rf so that 


and 


Hence, 


(a\u k \ 2 W^ * \auo k \ 2 ) (-,- 2 ) = - (alcjfcl 2 !!^ 7 * [awfcl 2 ) (-,z) 
f a(x)\u}k(z)\ 2 W 2 * \a(x)uk(z)\ 2 dz = 0 . 

J R d 

1 /* 27 r -1 r 2 n 

F 2M = —j^ F 1 (9,u)d0 + — J F^(9,u)d9. 


The first step of the proof of Lemma 2.3 gives us the result. 


□ 









DIMENSION REDUCTION FOR DIPOLAR. BEC 


15 


3. Proofs of our main Theorems 

This section is devoted to the proofs of Theorems 1.6 and 1.7 and Proposition 
1.11 which are inspired by the ones of [4, Theorem 1.3. and 1.4.]. To do so, we 
recall without any proof some of the results the authors obtained in this paper for 
the sake of readability. 

3.1. Main tools. We begin by the following Lipschitz estimates which summarize 
[4, Lemma 2.7.] and Lemma 2.2. 

Proposition 3.1. For all m >2, there exists C m > 0 such that 

II Fav( u ) - F 2v( v )\\b™ < C m M 2 \\u - v\\ B ™ 

\\F^(0,u)-F'v(9,v)\\ Brn < C m M 2 \\u — v\\ B ™, 

for all M > 0, u,v £ B m (R 3 ) satisfying ||rt||sm < M, ||u||s m < M, all 6 £ M and 
all 7 £ [0,1]. 

We give then in Proposition 3.2 the local in time well-posedness of the eikonal 
equation [4, Proposition 2.2.]. 

Proposition 3.2. If So £ SQ s+ i(M 3_rf ) with s > 2, there exists T > 0 such that 
the eikonal equation (1.18) admits a unique solution S £ C([0, T]\ SQ s (R 3-d )) H 
C s ([0, T] x R 3 ~ d ). 

The following lemma is related to the non-homogeneous linear equation (3.1) (see 
[4, Lemma 2.6.]). The crucial bound (3.2) is obtained by energy estimate. 

Lemma 3.3. Let us assume that for some m > 2, s>m + 2 and T > 0, we have 

(i) a 0 £ B m (R 3 ), 

(ii) S £ C([0,T]; SQ s (M 3 ~ d )) nC s ([0,T] x R 3-d ) solves the eikonal equation (1.18), 
(Hi) i?£C([0,T];S m (M 3 )). 

Then, for all a £ [0,1], there exists a unique solution a € C([0, T]; B m (R 3 )) n 
C 1 ([0, T]; B m ~ 2 (R 3 )) to the following equation: 

CL CX. 

d t a + V X S ■ V x a + -A X S = i—A x a + R, a(0, x, z) = a 0 (x, z). (3.1) 

Moreover for all t £ [0,T], a satisfies the estimates 

||a(t)||^ m < Haoll^m + C f ||a(s)||| m ds + f (a(s), R(s)) Bm ds (3.2) 

Jo Jo 

< \\ao\\ 2 B m + C J^ ^||a(s)||| m + ||i?(s)||sm^ ds (3.3) 

where C is a generic constant which depends only on m and on 

sup ||<9^|| i cx> ([0T]><K 3-d ) . 

2<N<s 

3.2. Proofs of Theorems 1.6 and 1.7 and Proposition 1.11. Theorem 1.6 can 
be proved by standard techniques. Point (i) is a consequence of Proposition 3.2. 
The existence and uniqueness result (ii) stems from a fixed-point technique based 
on the Duhamel formulation of the different equations and on the local Lipschitz 
estimates of Proposition 3.1. The uniform bound (iii) can be obtained by Gronwall 
lemma. For details, one can refer for instance to [4] where the case Ao = 0 was 
treated. 
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Let us now prove Theorem 1.7. 

Averaging limit £ —>• 0: proof of (1.27). For 7 G [0,1], let us introduce the function 

: K x B m (R 3 ) — > B m (R 3 ) 

( 0,u ) 1 —» Jq(F 1 (s,u) — Ff v {u))ds 

which satisfies the following properties for every u G B m (R 3 ): 

(a) 6 1 —^ is a 27r-periodic function, since 6 i->- F' y (d,u ) is 27r-periodic 

and Fa V is its average, 

(b) if ||rt||fjm < M then ||j^ 7 (0,u)\\b™ < 47 rC m M 3 for all 0 G R, where C m was 
defined in Proposition 3.1. 

Using the relations 

FHt/E 2 ,u(t)) = (F-<(t/E 2 Mt)) - F2M(t))) + n ,(«(*)), 

u(t))) = (f’WE 2 , »(*)) - e?„(u(t))) +E 1 D u ^(t/E l ,u(t))(d,u(t)), 

and equations (1.24) and (1.23), we obtain for all (a, 7 ) G [0, l ] 2 and e G (0,1], 

+ Va;5 ■ V x + ^ - y(A £ ’“’ 7 - A°’“’ 7 ) = (3.4) 

-i (F a 7 (A £ ’“’ 7 ) - F a 7 (A°’“’ 7 )) - i£ 2 d t ^(t/e 2 ,A £ ^) 

+ie 2 D u ^ (t/e 2 , A £ ’ a > 7 ) (<9 t A £ >“’ 7 ). 

We have that 

sup sup ||H n ^" 7 (s/e 2 , A £ ’ a, ' y (s))(dtA e,a, ' y (s))\\ B m -2 < C. (3.5) 

se[o,T] e - Q 

Indeed, according to Theorem 1.6, the sequences 


(A 


£,",71 


Je,a, 7 


and {d t A £ ’ a ^] 


e,a,7 


are uniformly bounded, respectively in 

L°°([0, T]; L> m (M 3 )) and in L°°([0, T]; B m - 2 (M 3 )). 


Thanks to Lemma 2.1, we get 



so that the dipolar part of Du^ 1 (t/e 2 , A e ’ a, ' y )(dtA £,a, ' Y ) is uniformly bounded in 
L°°([0, T]; i? m ^ 2 (M 3 )). The same property also holds for the cubic nonlinearity so 
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that we obtain inequality (3.5). Then, applying Lemma 3.3 to equation (3.4), we 
get 

\\A e ’ a " _ < (J f _ A°>“>l| m _ a (s)cfe 

JO 

+ f || FZ(A e ™) - F2M^)f Bm . 2 ds 

Jo 

+ e 4 [ ||D u ^(t/e 2 ,A £ '^)(9^ £ ’“^)||^ m _ 2 d S 
Jo 

+ e 2 f (d t ^(t/s 2 ,A £ ^),A £ ^ - A°>“") flm _ 2 (s)ds 

Jo 

<Ce A + C ( || A 6 ’*" - A°^\\ 2 Bm . 2 ( s)ds 

Jo 

- e 2 f {d t ^(t/e 2 ,A^),A^ - A°^) Bm _ 2 ( s)ds , 

Jo 

where we used the Lipschitz estimates of Proposition 3.1. The last term can be 
treated exactly as in [4], integrating by parts in time and using the equation (3.4). 
The conclusion follows by the Gronwall lemma. 

The semi-classical limit a —>• 0: proof of (1.28). The proof of the error estimate 

(1.28) follows exactly the same arguments as the ones of [4, Theorem 1.4.] since for 
7 fixed, the new dipolar term can be treated exactly as the cubic term. 

The dipole-dipole interaction limit 7 —>• 0: proof of (1.29). In the case of inequality 

(1.29) , we have for any e € (0,1], (a, 7 ) € [0, l ] 2 that 

d t (A £ ’ a ’7 - A £ ’ a ’°) + X7 x S ■ V x (A £ ’ a " - A £ ’ afi ) + =^(A £ ’ a ^ - A £ ^°) 

= i-A x (A £ ’ a - A £ ’ a ’°) - i(F' y (s/e 2 ,A £,a,J ) - F^(s/e 2 , A £ ’ Q ’ 0 )) 

-i (F 1 (s/e 2 , A £ ’ a, °) - F°{s/E 2 ,A £ ’ a ’ 0 )) 

so that Lemma 3.3, Proposition 3.1, Theorem 1.6 and Lemma 2.3 ensure that 


\A £ ’ a ^(t)-A‘ 

rt 


+ 


+ 


< 


Will™-. < c / - ^' 0 , '°MII|~. 

Jo 

f (II F^s/e 2 ,A £ ^(s)) - F^s/e 2 ,A £ ^°(s)) ds 

Jo 

f (II F^s/e 2 ,A £ ^°(s)) - F°(s/e 2 ,A £ ^°(s ))||| m _ 5 ) ds 

Jo 

C^ + C [ \\A £ ^(s) - A £ ’ a ’°(s)\\ 2 Brn - 5 ds 

Jo 


ds 


for any q € (0,1), and by Gronwall’s lemma 

11 /ie,a,7 _ /.°,a,7|| _ < (TV* 

II 71 71 IIC([0,T];B m - 5 ) — °7 ■ 

The case £ = 0 follows the same ideas. The proof of (1.29) is complete and (1.30) 
follows. □ 
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Proof of Proposition 1.11. In this case, we remark that the solutions remain po¬ 
larized on a single mode of PL Z as time evolves. Hence, we can apply Lemma 2.4 
instead of Lemma 2.3 and Proposition 1.11 follows from the arguments used in the 
proof of the estimate (1.29) in Theorem 1.7 
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